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PREFACE 


1.  The  Corps  of  Engineers,  through  its  Civil  Works  program,  lias  sponsored,  over  the  past 
23  years,  research  into  the  behavior  and  characteristics  of  tidal  inlets.  The  Corps’  interest  in 
tidal  inlet  research  stems  from  its  responsibilities  for  navigation,  beach  erosion  prevention 
and  control,  and  flood  control.  Tasked  with  the  creation  and  maintenance  of  navigable  l .S. 
waterways,  the  Corps  routinely  dredges  millions  of  cubic  yards  of  material  each  year  from 
tidal  inlets  that  connect  the  ocean  with  bays,  estuaries,  and  lagoons.  Design  and 
construction  of  navigation  improvements  to  existing  tidal  inlets  are  an  important  part  of  the 
work  of  many  Corps’  offices.  In  some  cases,  design  and  construction  of  new  inlets  are 
required.  Development  of  information  concerning  the  hydraulic  characteristics  of  inlets  is 
important  not  only  for  navigation  and  inlet  stability,  but  also  because  inlets,  by  allowing  for 
the  ingress  of  storm  surges  and  egress  of  flood  waters,  play  an  important  role  in  the  flushing 
of  bays  and  lagoons. 

2.  A research  program,  General  Investigation  of  Tidal  Inlets  (G1TI),  was  developed  to 
provide  quantitative  data  for  use  in  design  of  inlets  and  inlet  improvements.  It  is  designed  to 
meet  the  following  objectives: 

To  determine  the  effects  of  wave  action,  tidal  flow,  and  related  phenomena  on  inlet 
stability  and  on  the  hydraulic,  geometric,  and  sedimentary  characteristics  of  tidal 
inlets;  to  develop  the  knowledge  necessary  to  design  effective  navigation 
improvements,  new  inlets,  and  sand  transfer  systems  at  existing  tidal  inlets;  to  evaluate 
the  water  transfer  and  flushing  capability  of  tidal  inlets;  and  to  define  the  processes 
controlling  inlet  stability. 

3.  The  G1T1  is  divided  into  three  major  study  areas:  (a)  inlet  classification,  (b)  inlet 
hydraulics,  and  (c)  inlet  dynamics. 

a.  Inlet  Classification.  The  objectives  of  the  inlet  classification  study  are  to  classify 
inlets  according  to  their  geometry,  hydraulics,  and  stability,  and  to  determine  tire 
relationships  that  exist  among  the  geometric  and  dynamic  characteristics  and  the 
environmental  factors  that  control  these  characteristics.  The  classification  study  keeps  the 
general  investigation  closely  related  to  real  inlets  and  produces  an  important  inlet  data  base 
useful  in  documenting  the  characteristics  of  inlets. 

b.  Inlet  Hydraulics.  The  objectives  of  the  inlet  hydraulics  study  are  to  define  the 
tide-generated  flow  regime  and  water  level  fluctuations  in  the  vicinity  of  coastal  inlets  and 
to  develop  techniques  for  predicting  these  phenomena.  The  inlet  hydraulics  study  is  divided 
into  three  areas:  (1)  idealized  inlet  model  study,  (2)  evaluation  of  state-of-the-art  physical 
and  numerical  models,  and  (3)  prototype  inlet  hydraulics. 
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(1)  The  Idealized  Inlet  Model.  The  objectives  of  this  model  study  are  to  determine 
the  effect  of  inlet  configurations  and  structures  on  discharge,  head  loss  and  velocity 
distribution  for  a number  of  realistic  inlet  shapes  and  tide  conditions.  An  initial  set  of  tests 
in  a trapezoidal  inlet  was  conducted  between  1967  and  1970.  However,  in  order  that 
subsequent  inlet  models  are  more  representative  of  real  inlets,  a number  of  “idealized” 
models  representing  various  inlet  morphological  classes  are  being  developed  and  tested.  The 
effects  of  jetties  and  wave  action  on  the  hydraulics  are  included  in  the  study. 

(2)  Evaluation  of  State-of-the-Art  Modeling  Techniques.  The  objectives  of  this 
part  of  the  inlet  hydraulics  study  are  to  determine  the  usefulness  and  reliability  of  existing 
physical  and  numerical  modeling  techniques  in  predicting  the  hydraulic  characteristics  of 
inlet-bav  systems,  and  to  determine  whether  simple  tests,  performed  rapidly  and 
economicall) , are  useful  in  the  evaluation  of  proposed  inlet  improvements.  Masonboro  Inlet, 
North  Carolina,  was  selected  as  the  prototype  inlet  which  would  be  used  along  with 
hydraulic  and  numerical  models  in  the  evaluation  of  existing  techniques.  In  September  1969 
a complete  set  of  hydraulic  and  bathymetric  data  was  collected  at  Masonboro  Inlet. 
Construction  of  the  fixed-bed  physical  model  was  initiated  in  1969,  and  extensive  tests  have 
been  performed  since  then.  In  addition,  three  existing  numerical  models  were  applied  to 
predict  the  inlet’s  hydraulics.  Extensive  field  data  were  collected  at  Masonboro  Inlet  in 
August  1974  for  use  in  evaluating  the  capabilities  of  the  physical  and  numerical  models. 

(3)  Prototype  Inlet  Hydraulics.  Field  studies  at  a number  of  inlets  are  providing 
information  on  prototype  inlet-bay  tidal  hydraulic  relationships  and  the  effects  of  friction, 
waves,  tides,  and  inlet  morphology7  on  these  relationships. 

c.  Inlet  Dynamics.  The  basic  objective  of  the  inlet  dynamics  study  is  to  investigate 
the  interactions  of  tidal  flow,  inlet  configuration,  and  wave  action  at  tidal  inlets  as  a guide 
to  improvement  of  inlet  channels  and  nearby  shore  protection  works.  The  study  is 
subdivided  into  four  specific  areas:  (1)  model  materials  evaluation,  (2)  movable-bed 
modeling  evaluation,  (3)  reanalysis  of  a previous  inlet  model  study,  and  (4)  prototype  inlet 
studies. 

• 

(1)  Model  Materials  Evaluation.  This  evaluation  was  initiated  in  1969  to  provide 
data  on  the  response  of  movable-bed  model  materials  to  waves  and  flow  to  allow7  selection 
of  the  optimum  bed  materials  for  inlet  models. 

(2)  Movable-Bed  Model  Evaluation.  The  objective  of  this  study  is  to  evaluate  the 
state-of-the-art  of  modeling  techniques,  in  this  case  movable-bed  inlet  modeling.  Since,  in 
many  cases,  movable-bed  modeling  is  the  only  tool  available  for  predicting  the  response  of 
an  inlet  to  improvements,  the  capabilities  and  limitations  of  these  models  must  be 
established. 

I- 
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(3)  Reanalysis  of  an  Earlier  Inlet  Model  Study.  In  1975,  a report  entitled, 
“Preliminary  Report:  Laboratory  Study  of  the  Effect  of  an  Uncontrolled  Inlet  on  the 
Adjacent  Beaches,”  was  published  by  the  Beach  Erosion  Board  (now  CERC).  A reanalysis  of 
the  original  data  is  being  performed  to  aid  in  planning  of  additional  GITI  efforts. 

(4)  Prototype  Dynamics.  Field  and  office  studies  of  a number  of  inlets  are 
providing  information  on  the  effects  of  physical  forces  and  artificial  improvements  on  inlet 
morphology.  Of  particular  importance  are  studies  to  define  the  mechanisms  of  natural  sand 
bypassing  at  inlets,  the  response  of  inlet  navigation  channels  to  dredging  and  natural  forces, 
and  the  effects  of  inlets  on  adjacent  beaches. 

4.  This  study  supports  the  goals  of  the  General  Investigation  of  Tidal  Inlets  by  using  a 
movable-bed  inlet  model  to  study  inlet  hydraulics  for  a variety  of  inlet  configurations  and 
for  various  conditions,  including  wave  steepness  with  and  without  jetties,  freshwater  inflow, 
and  bed  rippling.  The  study  also  suggests  parameters  useful  to  classify  inlet  hydraulics  and 
briefly  discusses  inlet  stability  by  re-examining  the  inlet  cross-sectional  area  versus  prism 
relationsiiip. 
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Ke  . 

37  Dimensionless  timelag  between  low  waters  as  a function  of  Kg 

38  Dimensionless  maximum  inlet  velocity, as  a function  of  Kg.  . . . 

39  Dimensionless  bay  superelevation  as  a function  of  Kg  „ ...  . 

40  Typical  rippled  beds  after  ebbtide  flow  for  high  (F)  inlets  . . 

41  Examples  of  typical  rippled  beds  after  ebbtide  flow  for 

low  (F)  model  inlets  


42  Fixed  ripple  beds  and  estimated  friction  coefficients  at 
maximum  discharge.  ....  


43  Model  inlet  tidal  prism  versus  area  . 


CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 

UNITS  OF  MEASUREMENT 

U.S.  customary  units  of  measurement  used  in  tliis  report  can  be  converted  to  metric  (SI) 
units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

bee  lares 

acres 

0.4047 

hectares 

foot-pounds 

1 .3558 

newton  meters 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton, long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

'To  obtain  Celsius (G)  temperature  readings  from  Fahrenheit  (F)  readings,  use  formula:  C = (5/9)  (F  - 32). 
To  obtain  Kelvin  (K)  readings,  use  formula:  K = (5/9)  (F  - 32)  +273.15. 
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Ac 


Ace 

Ac/ 


A pitot 


ao 

ab 

C 

CD 

CDi 

Cf 

C9 

o 

D 

Dc 

d 
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F 

f 

g 

H 

HBL 

MIL 

HOL 


SYMBOLS  AND  DEFINITIONS 

bay  surface  area  (constant  for  vertical  walls)  in  square  feet 

inlet  channel  cross-sectional  area  at  any  time  (generally  at 
mean  sea  level,  MSL)  in  square  feet 

mean  channel  cross-sectional  area  (runs  A)  or  cross-sectional 
area  at  throat  (runs  B and  C)  below  MSL  at  maximum  ebbtide  and 
floodtide  discharges  (square  feet) 

pilot  channel  cross-sectional  area  below  MSL  (runs  A)  or  average 
area  of  preceding  run  below  MSL  for  relevant  runs  (runs  B and  C) , 
in  square  feet 

ocean  tide  amplitude,  equal  to  one-half  the  ocean  range 
bay  tide  amplitude,  equal  to  one-half  the  bay  range 
celerity  of  surface  waves  (feet  per  second) 
velocity  coefficient  at  maximum  discharge  (dimensionless) 
discharge  coefficient  of  maximum  discharge  (dimensionless) 

Chezy  friction  coefficient; (feet)'2  per  second 
wave  group  velocity  (feet  per  second) 
a subscript  meaning  channel 
water  depth  (feet) 

water  depth  below  MSL  near  the  current  meter 

hydraulic  depth  at  MSL  (feet) 

diameter  of  sand  grain  (millimeters) 

total  energy  in  one  wavelength  to  permit  crest  width 
(foot-pound  per  foot) 

inlet  energy  loss  coefficient  (dimensionless) 

Darcy-Weisbach  friction  coefficient  (dimensionless) 

acceleration  of  gravity  (feet  per  second  squared) 

average  wave  height  at  ocean  beach  toe  (feet) 

highest  bay,  inlet,  and  ocean  water  levels  referred  to  MSL 
(feet) 
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Ho 

h0 


h&e 
h bf 

hoe 
h of 


hoe 
h of 

h ie 
h if 

K 


hen 


SYMBOLS  AND  DEFINITIONS--Continued 
average  wave  height  in  deep  water  (feet) 

average  height  of  highest  ripples  along  channel  centerline 
(feet) 

instantaneous  water  surface  elevation  at  center  of  bay  (feet) 

instantaneous  water  surface  elpvation  in  ocean  at  toe  of 
beach  (feet) 

bay  elevations  at  maximum  ebbtide  and  floodtide  discharges 
(feet) 

water  surface  elevations  at  midchannel  (runs  A)  or  at  channel 
throat  (runs  B and  C)  at  maximum  ebbtide  and  floodtide 
discharges  (feet) 

ocean  elevation  at  maximum  ebbtide  and  floodtide  discharges 
(feet) 

inlet  water  level  at  maximum  ebbtide  and  floodtide  discharges 
(feet) 

Keulegan  repletion  coefficient,  dimensionless  (subscript  e 
for  ebbtide,  f for  floodtide) 

coefficient  of  energy  losses  due  to  channel  entrance 
(dimensionless) 


kSx  coefficient  of  energy  losses  due  to  channel  exit 

(dimensionless) 

L length  of  channel:  distance  along  centerline  from  start  of 

ripples  in  ocean  beach  to  crest  of  bay  delta;  used  for 
estimation  of  friction  coefficients  (feet) 

L0  average  wavelength  in  deep  water  (feet) 


MBL  mean  bay  elevation  above  MSL  over  a tidal  cycle  (feet) 


MIL  mean  inlet  water  level  at  throat  (runs  B and  C)  or  at 

midchannel  (runs  A)  over  a tidal  cycle  (feet) 


MSL  mean  ocean  level  over  a tidal  cycle.  The  horizontal  plane 

at  this  elevation  is  the  datum  to  which  all  vertical  distances 
are  referred.  The  value  of  MSL  in  Table  1 is  the  location  of  I 
MSL  below  the  Stillwater  level  (feet) 


Manning's  roughness  coefficient 
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SYMBOLS  AND  DEF INIT IONS- -Continued 


p average  wetted  perimeter  (runs  A)  or  wetted  perimeter  at  throat 

(runs  B and  C)  at  MSL  (feet) 

p0  average  wave  power  in  deep  water  (pound-foot  per  second  per 

foot  of  wave  crest) 

Q discharge  through  inlet  at  time  t (cubic  foot  per  second) 

"freshwater"  discharge  into  bay  (cubic  foot  per  second) 

Qmcuc  maximum  discharge  through  inlet  (cubic  foot  per  second) 

R hydraulic  radius  at  time  t (feet) 

R^  range  of  tide  at  center  of  bay  (feet) 

RMSL  hydraulic  radius  at  MSL  (feet) 

Rc  range  of  tide  at  toe  of  ocean  beach  (feet) 

S slope 

S be  rate  of  emptying  of  bay  at  maximum  ebbtide  discharge  and  of 

Sbf  filling  of  bay  at  maximum  floodtide  discharge  (0.00001  foot 

per  second) 

STAT  station  at  which  throat  is  located  according  to  reference  scale 
along  model  inlet;  equal  to  zero  for  runs  A 

SWL  Stillwater  level  common  to  ocean,  channel,  and  bay  at  the  end 

of  each  run 

T tidal  period:  duration  of  tidal  cycle  (seconds) 

TE  approximate  duration  of  ocean  drop  (seconds) 

T0  average  wave  period  (seconds) 

t elapsed  time  (seconds) 

V velocity  of  current  (feet  per  second) 

Vmax  measured  maximum  velocities  at  channel  centerline  (midchannel 

for  runs  A and  at  throat  for  runs  B and  C)  at  maximum  discharge 
(feet  per  second) 

^max  computed  maximum  average  velocity  across  the  average  cross 
section  (runs  A)  or  across  the  throat  (runs  B and  C)  (foot 
per  second) 

I dimensionless  mean  maximum  velocity  (equation  5)  across  the 

average  channel  cross  section  (rus  A)  or  across  the  throat 


SYMBOLS  AND  DEFINITIONS--Continued 


W MSL 
w 

max  or  m 
min 
b 
e 

f 

T 

o 

a. 

e 

y 
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average  channel  width  (runs  A)  or  channel  width  at  throat 
(runs  B and  C)  at  MSL  (feet) 

width  of  inlet  channel  at  a given  location  at  time  t (feet) 

subscript;  maximum 

subscript;  minimum 

subscript;  bay  condition 

subscript;  at  maximum  ebbtide  flow 

subscript;  at  maximum  floodtide  flow 

superscript;  theoretical 

subscript;  deepwater  condition  (ocean)  or  MSL  condition 

velocity  distribution  coefficient  (dimensionless) 

coefficient  for  linear  variation  of  horizontal  bay  surface 
area  (dimensionless) 

unit  weight  of  water  (pounds  per  cubic  foot) 

superelevation  of  bay  over  ocean  defined  as  the  difference  in 
elevations  between  MBL  and  MSL  over  a tidal  cycle  (feet) 


£ average  channel  bank  slope  (runs  A)  or  bank  slope  at  throat 

(runs  B and  C) ; horizontal  distance  per  unit  of  vertical 
raise 

dimensionless  timelags  between  ocean  and  bay  high  or  low 
waters  = timelag  per  tidal  period 

n water  surface  elevation  at  inlet  taken  as  the  average  of 

instantaneous  bay  and  ocean  levels  (feet) 

Q tidal  prism  equal  to  R ^ times  A^  for  all  runs  (cubic  feet) 
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LABORATORY  INVESTIGATION  OF  TIDAL  INLETS  ON  SANDY  COASTS 

by 

Ramiro  E.  Mayor-l-iora 
I.  INTRODUCTION 

Sandy  beaches  are  a common  feature  along  the  coast  and  often  act  as 
barriers  between  the  ocean  and  inland  bodies  of  water  (i.e.,  bays,  la- 
goons, estuaries).  Such  barriers  may  have  one  or  more  relatively  narrow 
channels  (inlets)  that  connect  the  two  bodies  of  water  and  control  the 
amount  and  the  manner  in  which  the  energy  of  ocean  tides,  waves,  and 
currents  is  exchanged  within  the  ocean-inlet-bay  system. 

Meteorological,  hydrographical,  geographical,  and  wave  conditions 
may  force  inlets  to  migrate,  change  their  orientation  or  shape,  and  at 
times,  close  completely.  Man's  intervention  in  the  natural  process  of 
an  inlet,  bay,  or  on  the  neighboring  coast  may  affect  an  inlet's  migra- 
tional  tendency  to  close. 

A study  of  the  hydrodynamic  phenomena  involved  is  complicated  by  the 
many  factors  that  act  upon  the  system  and  by  the  unpredictable  nature  of 
most  of  these  factors.  This  study  is  designed  to  provide  a better  under- 
standing of  the  factors  influencing  inlet  behavior  which  may  affect  nav- 
igation, water  quality,  the  ecology,  recreation,  and  marine  life.  Both 
the  dynamic  changes  in  the  inlet  bed  and  the  nearby  beaches  and  the 
hydraulic  characteristics  of  the  system  are  major  aspects  of  the  problem 
that  many  investigators  have  studied  in  recent  years. 

Si  the  collection  of  field  data  presents  many  problems,  completely 
controlled  model  studies  are  needed  to  examine  basic  inlet  processes  and 
parameters.  In  response  to  this  need  for  additional  understanding  of 
general  inlet  processes,  this  study  includes: 

(a)  Laboratory  testing  of  a movable-bed  inlet  with  varying 
geometric  characteristics  subjected  to  combinations  of  ocean 
tides  and  waves; 

(b ) measurement  of  the  temporal  and  spatial  variation  in 
water  levels,  current  velocities,  and  channel  geometry; 

(c)  investigation  of  the  effects  of  tide,  wave,  runoff, 
and  geometric  parameters  on  the  inlet  channel  hydrography  and 
the  flushing  ability  of  the  system; 

(d)  correlation  between  dimensionless  hydraulic  parameters; 

and 

(e)  comparison  of  model  results  to  theoretical  hydraulic 
solutions . 
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This  laboratory  study  uses  several  simplifications  to  reduce  the 
complexity  of  the  inlet  characteristics,  processes,  and  data  analysis. 

These  include:  A fixed  inlet  location;  no  study  of  the  stability  of 

the  channel  nor  the  sediment  transport  rates  through  it;  no  consideration 
of  saltwater  effects  or  flow  patterns  at  the  entrance  and  vicinity;  and 
no  investigation  of  longshore  transport. 

The  report  presents  a brief  review  of  the  principal  theoretical 
approaches  developed  for  predicting  the  hydraulic  characteristics  of  the 
system,  a summary  of  field  and  other  laboratory  data,  and  a description 
of  the  experimental  program.  The  collected  data  are  given  in  reduced 
form  and  the  various  dimensionless  parameters  discussed,  correlated,  and 
compared  to  theoretical  approximations.  Comparisons  of  laboratory  and 
field  data  are  included.  Conclusions,  comments,  and  suggestions  for 
further  use  of  the  empirical  data  are  also  presented. 

II.  THEORETICAL  ANALYSES 

i 

1 . Simplified  Tidal  Inlet  Hydraulics. 

A complete  study  of  the  hydraulics  of  a bay- inlet- ocean  system  would 
require  consideration  of  the  numerous  time- dependent  factors  acting 
upon  such  a system,  including: 

(a)  Ocean  climate  (wind  waves,  storm  surge,  tide); 

(b)  nearshore  conditions  (longshore  currents,  manmade 
modifications  to  the  shore); 

(c)  hydrologic  conditions  in  the  basins  draining  into 
the  bay; 

(d)  geometric  characteristics  of  the  bay  (shape,  depth, 
irregularity,  slope  of  sides)  and  the  inlet  channel  (location 
relative  to  the  bay,  depth,  orientation,  length,  shape,  and 
ability  to  shift  location); 

(e)  inlet  sediment  characteristics;  and 

(f)  ocean- freshwater  density  wedge  effects. 

Since  these  factors  are  beyond  the  scope  of  this  investigation,  the 
study  has  been  simplified  to  the  determination  of  water  level  fluctuations 
in  the  bay  and  flow  through  the  inlet  channel  (average  current  velocity, 
discharge)  in  response  to  a defined  periodic  ocean  tide.  Therefore,  the 
approaches  found  in  the  literature  and  discussed  in  this  section  involve 
some  simplifying  assumptions  common  to  all  of  these  analyses:  relatively 

small  and  deep  bays;  constant  water  density  throughout  the  system;  fixed 
sediment  characteristics  and  behavior;  and  a fixed  inlet  channel  orien- 
tation, location,  and  shape. 
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2.  Gradually  Varied  Unsteady  Flow  in  Open  Channels. 


The  simplified  equations  governing  gradually  varied  unsteady  flow  of 
water  in  an  open  prismatic  channel  are  the  continuity  equation: 


D 


3V 

3x 


+ 


v 1Z.  + *L  o o 
3x  3t 


(1) 


and  the  general  dynamic  equation: 


3y  V 3V  1 3V 

3x  g 3x  g 3t 


S(?  + S = 0 


(2) 


where  x and  y define,  respectively,  a horizontal  axis  along  the  channel 
and  a vertical  axis  along  which  heights  above  a datum  are  measured;  D is 
the  hydraulic  depth  (A<?/w) , where  Ac  is  the  cross-sectional  area  and 
w is  channel  width;  g is  the  acceleration  due  to  gravity;  3y/3x  and 
3y/3t,  the  variation  of  the  water  surface  with  distance  and  with  time, 
respectively;  Sc  is  the  slope  of  the  channel  bottom  along  its  longi- 
tudinal axis;  S is  the  slope  of  the  energy  line;  and  a is  the  coef- 
ficient resulting  from  the  velocity  distribution  (Chow,  1959). 

The  change  in  total  head  for  the  channel  depends  on  friction,  accel- 
eration, and  eddy  (turbulence)  losses. 

3.  Principal  Theoretical  Approaches  to  the  Problem. 

The  application  of  equations  (1)  and  (2)  to  an  ocean-inlet-bay 
system  has  occupied  the  attention  of  several  investigators  (Chapman, 
1923;  Brown,  1928;  O'Brien,  1937;  Keulegan,  1967;  Van  de  Kreeke,  1967; 
Shemdin  and  Forney,  1970;  Mota  Oliveira,  1971;  Huval  and  Wintergerst, 
in  preparation,  1977).  The  main  differences  among  these  studies  are 
the  assumptions  regarding  the  physical  phenomena  and  the  relative 
importance  given  to  the  basic  terms  included  in  these  equations.  The 
principal  theoretical  approaches  are  summarized  below. 

a.  The  Keulegan  Solution.  Keulegan  (1967)  obtained  a mathematical 
method  of  describing  inlet  hydraulics  with  these  assumptions: 

(a)  Walls  of  the  bay  are  vertical; 

(b)  there  is  no  inflow  from  streams; 

(c)  no  density  currents  are  present; 

(d)  tidal  fluctuations  are  sinusoidal; 

(e)  the  bay  water  level  rises  uniformly; 
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(£)  length  of  the  bay  is  less  than  the  tidal  wavelength; 

(g)  the  inlet  channel  flow  area  is  constant;  and 

(h)  the  inertia  of  the  mass  of  the  water  in  the  channel 
is  negligible. 

An  ocean-inlet-bay  system  which  satisfies  the  Keulegan  assumptions  is 
shown  in  Figures  1 and  2. 

The  Keulegan  solution  describes  the  water  surface  fluctuation  in  the 
bay  by  the  differential  equation: 


dhb 

dt 


K /2  1 hD  - h^j  | » 


(3) 


where  h^  is  the  water  level  in  the  bay;  h0  is  the  water  level  of  the 
ocean;  T is  the  tidal  period;  R0  is  the  tidal  range  in  the  ocean;  and 
K is  the  coefficient  of  repletion,  a dimensionless  number  defined  as: 


K 


T Ae  I g-~  > 

tt/r ^ A£>  \ ken  + kex  + fL/4R 
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where  Ag  is  the  inlet  cross-sectional  area;  A^,  is  the  area  of  the  bay; 
kgn  is  the  entrance  loss  coefficient;  kgx  is  the  exit  loss  coefficient; 
f is  the  Darcy-Weisbach  friction  coefficient  of  the  channel;  L is  the 
channel  length;  and  R is  the  hydraulic  radius  of  the  channel. 

O'Brien  and  Dean  (1972)  presented  Keulegan's  solution  in  dimension- 
less bay  tidal  range  and  phase  lag  plotted  against  the  coefficient  of 
repletion,  K (Fig.  3).  Note  in  the  figure  that  for  low  values  of  K, 
the  tide  amplitude  in  the  bay,  a£,  is  smaller  than  the  ocean  amplitude, 
aa  (a£/a0  is  less  than  1).  As  K increases,  the  amplitude  in  the  bay 
increases,  so  that  for  values  of  K greater  than  2,  the  ocean  and  bay 
tide  amplitudes  are  the  same. 

Similarly,  the  lag  between  the  ocean  and  bay  water  levels  e, 
•approaches  90°  for  small  values  of  K and  drops  to  zero  for  values 
of  K greater  than  10  (Fig.  3). 

Dimensionless  maximum  velocity,  V,^.,  is  also  given  as  a func- 
tion of  K with  maximum  velocity  increasing  for  higher  values  of  K 
(Fig.  3).  Maximum  inlet  velocity,  Vrnax>  is  related  to  dimensionless 
maximum  velocity  by  the  equation: 


, a AJ 

^ max  ~ ^ max  T °A c 
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b Lumped  Parameter  Approach.  The  Keulegan  (1967)  method  of  pre- 
dicting inlet  hydraulics  makes  several  restrictive  assumptions,  i.e.,  a 
sinusoidal  tide,  vertical  bay  walls,  and  vertical  channel  walls.  The 
lumped  parameter  approach  (Huval  and  Wintergerst,  in  preparation,  1977) 
generalizes  the  Keulegan  solution  by  allowing  for  a wider  range  of 
conditions : 

(a)  A nonsinusoidal  ocean  tide; 

(b)  variable  bay  surface  area; 

(c)  variable  inlet  depth  and  cross-sectional  area 
throughout  a tidal  cycle  due  to  sloping  inlet  banks  and 
variable  water  levels; 

(d)  freshwater  inflow  into  the  bay;  and 

(e)  inclusion  of  the  acceleration  terms  of  equation  (2) . 

In  equation  (2),  the  dynamic  equation  for  flow  in  a channel  with  a 
bed  slope  of  zero  (So  = 0)  connecting  an  ocean  and  a bay  can  be 
simplified  to: 
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Taking 

Jay  = hb  - h0,  Jsx  = L,  j V3V  = ~ , 
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(quadratic  resistance  assumption)  and  rearranging, 


ho  - hfe  = — \ken  * kex  + |v|v  + | f^r  > 
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(7) 


where  h^  is  any  function  of  time. 

Introducing  the  continuity  equation  with  freshwater  inflow, 
(any  function  of  time) : 
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at  Ah 


24 


wn 


Equations  (7)  and  (8) 
varies  approximately  with 
the  term  | V | V,  both  Aj 


are  nonlinear.  Since  the  hydraulic  radius,  R, 
the  water  level  R - R^^  + p,  and  because  of 
and  are  also  functions  of  time.  Let: 


and 


N?  = A i)o  1 
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where  n = (h£  + h0)/2,  Afro  is  the  bay  surface  area  at  mean  sea  level 
(MSL) , g is  a dimensionless  parameter  of  bay  slope,  and  t,  is  the 
channel  wall  slope.  Furthermore,  the  friction  factor,  f,  is  expressed 
by  the  Manning  coefficient,  n,  and  the  hydraulic  radius  (Chow,  1959)  as: 


f = JtiLiL2  _ 1 16 . 7n2  (1]) 
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A Runge-Kutta-Gill  numerical  method  for  integration  of  the  two  simulta- 
neous ordinary  differential  equations  (7)  and  (8)  was  programed  for  the 
digital  computer,  taking  ken  + kgx  = 1.  Starting  with  initial  boundary 
conditions,  the  program  integrates  the  equations  over  several  tidal  cycles 
until  bay  tide  and  velocity  become  periodic. 

The  original  computer  program  (Huval  and  Wintergerst,  1977)  computes 
instantaneous  bay  tide,  average  inlet  current  velocity,  and  inlet  dis- 
charge for  given  values  of  Manning's  n.  This  program  was  modified  so 
that  the  expression  (ken  + kex  + fL/4R)  in  equation  (7)  is  substituted 
by  F,  a factor  determined  from  experimental  results,  thus  avoiding  the 
necessity  of  assuming  values  of  ken,  kex,  and  f. 

III.  THE  EXPERIMENTAL  PROGRAM 

1 . Experimental  Conditions. 

A systematic  series  of  36  laboratory  experiments  was  carried  out  on 
an  idealized  ocean-inlet-bay  system  involving  reversing  flows  (tidal 
action)  and  wave  action.  Del  Monte  white  sand,  60  mesh  with  a median 
diameter  of  0.34  millimeter,  was  used  to  model  the  barrier  in  all  tests 
(Fig.  4).  Freshwater  was  used. 

Other  parameters  constant  throughout  the  tests  were: 

(a)  The  surface  area  of  the  bay  (vertical  tank  walls) ; 

(b)  the  water  depth  at  MSL  in  the  ocean  basin; 

(c)  the  sand  barrier  slope  on  the  bay  side  and  the  bank 

slopes  of  the  pilot  channels; 
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Figure  4.  Sieve  analysis  of  the  sand  used  in  the  testi 


(d)  the  location  and  orientation  of  the  wave  generator; 

(e)  the  ocean  beach  slope  (1:30); 


(f)  the  alinement  of  the  pilot  channels  (parallel  to  the 

direction  of  wave  propagation) ; and 

(g)  the  relative  location  of  the  walls  defining  the 

ocean  basin,  the  sand  barrier  section,  and  the  bay  basin. 

The  independent  variables  (ocean  tidal  ijange,  pilot  channel  area  and 
length,  wave  characteristics,  tidal  cycle  amplitude  and  duration)  were 
selected  so  that  the  range  of  repletion  coefficients,  K,  varied  from 
0.04  to  1.30  to  provide  a wide  range  of  a^/a 0 ratios  (0.025  to  0.95). 
Experimental  conditions  tested  are  shown  in  Tables  1 and  2. 

In  addition  to  a series  of  tests  with  an  unimproved  inlet,  two  runs 
were  conducted  to  study  the  response  of  an  inlet  improved  by  parallel 
vertical-wall  jetties.  Further  studies  included  a run  with  "freshwater" 
discharge  into  the  bay,  runs  with  a "fixed-bed"  prismatic  channel,  and  a 
run  with  the  bed  form  "fixed"  in  plaster. 

This  study  does  not  include  a systematic  account  of  water  circulation 
patterns  in  the  ocean,  channel  or  bay,  or  any  measurement  of  sediment 
transport  rates  or  size  composition  of  the  bed;  it  concentrates  more  on 
measurements  of  water  surface  elevations  and  current  velocities  as  a 
function  of  time,  timelags,  and  the  size  and  shape  of  the  inlet  at  the 
end  of  each  run.  Temperature  effects  were  not  investigated,  although  at 
times  the  water  temperatures  in  the  ocean  and  in  the  bay  differed  by  8° 
Fahrenheit. 

2.  Testing  Facilities  and  Equipment. 

The  experiments  were  conducted  at  the  Richmond  Engineering  Field 
Station,  University  of  California,  Berkeley.  A 2.5-foot-deep,  94.3- 
foot-  long,  and  63.7-foot-wide  basin  was  divided  into  three  subbasins  by 
concrete  block  walls  covered  with  plastic  sheets  on  both  sides  to  prevent 
any  leakage  through  the  walls  (Fig.  5).  The  larger  subbasin  represented 
the  bay  and  was  separated  from  the  smaller  basin  (the  ocean)  by  a sand 
barrier  contained  between  two  guide  walls  (Fig.  6).  This  8- foot-wide 
sand  barrier  extended  into  the  ocean  with  a slope  of  1:30;  on  either 
side,  8-foot-wide  cement  barriers  formed  the  ocean  beach.  The  third 
basin  collected  the  water  outflow  from  the  ocean  where  it  drained  into 
the  outside  storage  basin. 

The  storage  basin,  located  outside  of  the  building,  was  equipped 
with  a filter  and  a pump.  The  circulating  system  operated  as  shown  in 
Figure  5. 

A 24- foot  flap- type  wave  generator  was  always  at  the  position  shown 
in  Figures  5 and  6,  so  that  the  periodic  waves  traveled  with  their  crests 
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perpendicular  to  the  beach.  Two  parallel  18-foot-long  steel  barriers, 
raised  about  6 inches  above  the  basin  bottom  to  allow  water  circulation, 
helped  prevent  spreading  of  wave  energy  caused  by  diffraction  and  refrac- 
tion. Wave  absorbers  were  installed  along  the  outer  walls  and  behind  the 
wave  generator. 

A special  mechanism  for  the  generation  of  sinusoidal  tides  was  de- 
signed and  constructed  (Fig.  7).  An  adjustable  10-foot-long  wooden  weir, 
with  a steel  angle  top  and  rubber  ends,  was  supported  by  three  turnbuckles 
connected  to  a common  shaft.  This  shaft  was  connected  to  a rotating  wheel 
through  a wooden  Pitman  arm  and  driven  by  a motor  equipped  with  a gearbox 
that  allowed  control  of  the  tidal  periods  by  varying  the  number  of  rev- 
olutions per  minute  (RPM).  The  driving  arm  could  be  fixed  to  the  wheel 
at  several  eccentric  positions  according  to  the  tidal  range  desired 
(about  equal  to  the  vertical  travel  of  the  weir).  The  turnbuckles  al- 
lowed further  adjustment  of  the  tidal  range,  the  mean  water  depth,  and 
leveling  of  the  weir  crest. 

The  inlet  region  and  the  moving  carriage  supported  by  30- foot- long 
rails  are  shown  in  Figure  8.  Instruments  mounted  on  the  carriage  reached 
any  point  on  the  sand  barrier.  Points  within  this  area  were  located  by 
perpendicular  scales  (in  feet)  mounted  on  a carriage  and  on  one  of  the 
steel  channels  supporting  the  rails.  The  scale  on  the  steel  channel 
defined  the  location  of  the  stations  at  which  measurements  were  made. 
Station  10  was  located  at  the  ocean  face  of  the  wall  dividing  the  two 
main  basins  and  station  0 was  on  the  bay  side.  A bridge  over  the  width 
of  the  basin  permitted  overhead  photography  (Fig.  8). 


Simultaneous  recording  of  water  surface  elevations  in  the  ocean  and 
current  velocities  in  the  channel  aided  control  of  the  tests.  As  shown 
in  Figure  9,  ocean  and  bay  tides  were  recorded  by  two  systems:  (a) 
Resistance- type  water  level  gages  (Fig.  10)  connected  to  two  Sanborn 
four-channel  recording  units  (Fig.  11);  and  (b)  Stevens  float-type 
water  level  recorders,  each  with  its  own  chart  mounted  on  a rotating 
drum  (Fig.  12).  Each  gage  and  each  float  was  protected  from  surface 
disturbances  by  appropriate  stilling  wells  with  openings  near  or  on  the 
basin  floor. 

Waves  were  measured  by  the  resistance- type  gages  and  recorded  by  a 
separate  two- channel  unit  with  varying  chart  speeds. 

A special  resistance  gage  was  used  to  detect  water  levels  at  the 
inlet  channel  because  of  its  relative  shallowness.  The  L-shaped  gage 
(mounted  on  the  carriage  as  shown  in  Fig.  13)  was  made  of  thin  wires 
spaced  0.25  inch  apart.  Three  other  special  gages  were  installed  along 
the  channel  axis  to  measure  the  water  surface  slope.  A standard  point 
gage  and  a small  current  meter  (Kent  Miniflow,  265-3,  propeller  diameter 
0.95  centimeter;  Fig.  14)  were  also  mounted  on  the  carriage.  The  current 
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Instrument  carriage  over  sand  barrier  and  platform  for  overhead 
photography. 
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Figure  9.  Schematic  layout  of  testing  facilities 

showing  relative  locations  of  instruments, 
controls,  and  recording  units. 
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meter  and  the  shallow-water  gages  were  connected  to  the  four- channel  re- 
corder units  to  record  ocean  tide,  bay  tide,  inlet  tide,  and  current 
velocity  on  the  same  chart. 
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Two  micromanometers  (readings  to  0.02  millimeter)  were  installed 
close  to  the  float  recorders  in  the  ocean  and  at  the  center  of  the  bay. 
Readings  at  high  and  low  tides  were  made  during  the  last  cycles  of  a run 
and  the  ranges  compared  with  those  from  the  float  charts.  These  readings 
proved  invaluable  when  water  level  data  were  reduced. 

The  two  temporary  float  recorders  (Fig.  9)  were  used  in  conjunction 
with  the  recorders  at  the  center  and  the  back  of  the  bay  to  check  the 
propagation  of  the  tidal  wave  in  the  bay. 

The  method  of  cross-sectioning  the  inlet  channel  after  the  end  of  a 
run  is  illustrated  in  Figure  15.  A vertical  planigraph  mounted  on  the 
bay  side  of  the  carriage  was  used  to  reproduce  the  sand  surface  elevation 
on  tracing  paper. 

4.  Testing  Program  and  Procedures. 

The  main  objective  of  the  tests  was  to  measure  the  hydraulic  response 
of  the  bay  and  movable-bed  inlet  to  tides  and  waves  in  the  ocean,  and  to 
the  initial  inlet  conditions. 

Since  the  development  of  equilibrium  in  a tidal  inlet  model  requires 
a long  testing  time  several  parameters  were  kept  constant  to  reduce  the 
number  of  tests.  The  constant  parameters  were:  (a)  The  slope  of  the 

ocean  beach  (1:30);  (b)  the  beach  material  (d$Q  = 0.34-millimeter  sand); 

(c)  the  geometry  and  surface  area  of  the  bay;  and  (d)  the  water  depth 
in  the  ocean  below  MSL  (0.91  to  0.94  foot). 

a . Character  of  Studies  and  Range  of  Independent  Variables.  Mos t 
tests  were  conducted  on  the  uncontrolled  (natural)  condition  of  the 
beach  and  inlet,  first  run  with  tidal  action,  then  re-run  with  simul- 
taneous tidal  and  wave  actions. 

Other  runs  provided  data  for  the  following  studies: 

(a)  Inlet  under  controlled  conditions  (jetties); 

(b)  effect  of  duration  of  run; 

(c)  propagation  of  the  ocean  tide  in  the  bay  (simultaneous 

vertical  movement  of  water  surface  at  different  locations); 

(d)  effect  of  initial  conditions  (pilot  channel  area 

and  water  level  in  bay) ; 

(e)  effect  of  wave  steepness; 

(f)  development  and  movement  of  bed  ripples; 
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(g)  effect  of  a fixed  prismatic  channel;  and 

(h)  effects  of  fixed  rippled  beds. 

The  channel  centerline  profile  and  inlet  cross-sectional  area  were 
measured  after  the  end  of  the  last  ebbtide  flow  and,  in  a few,  after 
floodtide. 

The  independent  variables  with  selected  ranges  were: 

(a)  Period  of  ocean  tidal  cycle*  (20  to  60  minutes), 

(b)  range  of  ocean  tide  (0.054  to  0.150  foot), 

(c)  pilot  channel  area  (0.082  to  0.434  square  foot), 

(d)  length  of  inlet  channel  (7  to  22.5  feet), 

(e)  wave  period  (0.75  to  0.95  second), 

(f)  wave  height  at  beach  toe  (0.030  to  0.176  feet), 

(g)  wave  steepness  (0.011  to  0.064), 

(h)  duration  of  run  (4  to  25  hours), 

(i)  water  temperature  (45°  to  76°  Fahrenheit)  uncontrolled, 

(j)  freshwater  inflow  discharge  (0.0  to  10.2  gallons 
per  minute) , and 

(k)  slope  of  bay  side  of  barrier  (1:3  to  1:12). 

b.  Preliminary  Runs.  Runs  2 and  3 were  devoted  to  preliminary  checks 
of  the  facilities,  i.e.,  testing  procedure,  equipment,  and  instruments. 

In  these  runs  the  inlet  barrier  was  formed  by  loosely  placing  just  enough 
dry  sand  to  separate  the  two  basins  (run  2)  and  form  a 7-foot-wide  barrier 
(run  3) . The  runs  helped  define  the  magnitude  of  the  independent  vari- 
ables used  in  subsequent  runs. 

c.  Preparations  for  a Typical  Run.  The  first  step  was  the  formation 
of  the  ocean  beach,  the  sand  barrier  crest,  and  the  bay-side  beach.  The 
basins  were  filled  with  freshwater,  about  2 inches  below  the  barrier 
crest,  to  saturate  the  sand.  A wooden  straightedge  supported  on  the 
interior  cement  beach  walls  was  then  progressively  moved  to  obtain  a 
flat  beach  with  a 1:30  slope.  The  same  procedure  was  followed  for  the 
crest  and  the  bay  beach,  taking  care  not  to  cause  any  compaction  of  the 
material.  Hxcess  sand  was  subsequently  removed. 

Cutting  of  the  pilot  channel  followed.  A trapezoid-shaped  board 
was  attached  to  the  front  of  the  carriage  and  pushed  along  the  inlet’s 
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longitudinal  axis  to  obtain  a centered  horizontal  prismatic  channel  with 
a base  1.70  feet  wide  and  30°  bank  slopes.  This  cutting  was  done  slowly 
to  avoid  accidental  wave  generation  and  resulting  bank  erosion.  Excess 
sand  was  removed  during  the  channel  cutting. 

Next  a uniform  Stillwater  level  (SIVL)  was  established  over  the  entire 
basin  using  care  to  avoid  flow  in  the  inlet. 

All  recording  units  were  balanced,  the  instruments  calibrated,  and 
the  SIVL  recorded  on  charts  and  floats  for  a short  time.  The  SWL  was 
also  read  from  micromanometers  and  point  gages.  The  instrument  carriage 
was  fixed  at  midchannel  before  the  actual  start  of  the  run. 

Three  series  of  test  conditions  used  during  the  investigation  were 
(prefix  n is  the  run  number) : 

Series  A (run  nA) ; tidal  action  only. 

Series  B (run  nB) ; tidal  and  wave  action. 

Series  C (run  nC) ; tidal  and  wave  action  with  special 

conditions . 

d.  Tidal  Action  Only  (Runs  nA).  The  resistance  gage  recording 
units  and  the  float  recorders  were  syncronized  (before  start  of  run)  by 
turning  on  the  tide  generator.  Normally,  the  weir  was  initially  at  the 
mean  position  and  started  its  upswing  movement  thus  creating  a floodtide 
current  into  the  bay.  Where  a superelevation  in  ^the  bay  was  expected 
(depending  on  barrier  and  inlet  characteristics,  tidal  cycle  and  tidal 
range),  the  bay  water  level  was  maintained  at  an  elevation  higher  than 
the  ocean  before  the  start  of  the  run.  This  was  done  by  temporarily 
placing  a thin  metal  blade  (Fig.  8)  across  the  channel  which  helped  to 
prevent  exaggerated  erosion  in  the  channel  during  the  first  floodtide 
flow  into  the  bay. 

High  rates  of  sediment  transport  were  visually  observed  during  the 
first  cycles,  most  of  it  as  bedload.  Ripples  were  formed  and  reversed 
orientation  and  movement  according  to  the  current  direction. 

Recording  of  water  stages  was  continuous  throughout  the  run.  Slack 
waters  were  marked  on  the  charts  by  manually  operating  a remote-control 
marker  when  reversal  of  the  flow  in  the  channel  was  observed.  The 
manometer  deflection  at  the  orifice  meter  in  the  inflow  pipe  was  con- 
stantly checked  to  ensure  that  the  discharge  was  constant  throughout  the 
run. 

During  the  last  tidal  cycles,  readings  of  high  and  low  water  eleva- 
tions were  taken  from  the  micromanometers,  current  velocities  were  re- 
corded, and  the  temperatures  of  the  bay  and  ocean  waters  were  recorded 
at  1-hour  intervals. 
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After  periodic  bay  tide  motion  was  attained,  the  metal  blade  was 
dropped  vertically  on  the  sand  inlet  at  the  slack  following  the  ebbtide 
flow,  to  isolate  the  bay  from  the  tidal  action.  The  tide  generator 
operated  until  the  weir  returned  to  the  initial  position;  then  the  circu- 
lation was  from  ocean  to  storage  basin  back  to  the  ocean.  The  pipe  be- 
tween the  bay  and  ocean  was  opened  to  obtain  a common  water  level  in  the 
two  basins  to  prevent  any  further  flow  through  the  inlet. 

All  recorders  continued  to  register  the  new  common  level  (,final  SKL, 
below  the  initial  level  in  all  runs)  obtained  when  the  blade  was  removed 
and  no  flow  was  observed  through  the  inlet  channel.  The  resistance  gages 
were  recalibrated  and  the  SWL  read  from  the  micrcmanometers  and  the  point 
gages. 

Finally,  the  pump  was  stopped  and  the  ocean  and  bay  drains  opened  to 
allow  the  slow  drainage  of  the  basins  into  the  storage  basin  until  a 
water  level  a few  inches  below  the  deepest  point  in  the  channel  bed  was 
attained.  In  this  procedure,  the  sand  bed  was  wet  at  all  times  and  the 
sand  grains  at  the  surface  layers  were  held  together  by  surface  tension 
and  capillary  forces,  an  ideal  condition  for  the  reproduction  of  the 
cross  sections  by  the  planigraph  (Fig.  15).  Normally,  this  cross- 
sectioning  was  made  at  stations  spaced  at  0.5  foot  on  the  day  following 
the  run. 

The  last  step  was  to  profile  the  channel  bottom  along  its  center- 
line  extending  from  the  bay  delta  to  a point  oceanward  of  the  original 
beach  unmodified  by  the  currents.  Vertical  photos  of  the  inlet  channel 
area  were  taken  from  the  overhead  platform. 

Durations  of  runs  A ranged  from  4 to  7.5  hours,  depending  on  the 
time  required  to  obtain  periodic  motion  of  the  bay  tides. 

Time-lapse  photos  in  Figure  16  illustrate  the  channel  bed  develop- 
ment (movement  of  sand  ripples)  for  a long  inlet  subjected  to  tidal 
action  only. 

e.  Simultaneous  Tidal  and  Wave  Action  (Runs  B and  C) . Runs  B were 
carried  out  after  part  A measurements  were  completed,  the  rippled  chan- 
nel becoming  the  initial  condition  for  the  new  run.  The  magnitude  of 
the  other  independent  variables  in  A remained  unchanged;  the  only  differ- 
ence was  the  introduction  of  periodic  waves. 

Filling  of  the  basins  and  preparations  for  the  run  were  similar  to 
those  described  previously.  If  bay  superelevations  were  detected  in 
previous  runs,  the  bay  level  was  initially  set  at  that  elevation  and  the 
metal  blade  removed  only  when  the  level  in  the  ocean  was  high  enough  to 

initiate  flooding  of  the  bay.  The  run  was  then  started,  and  the  selected 

waves  were  introduced  at  mean  ocean  level  only  after  one  or  two  cycles  of 

tidal  action.  This  procedure  was  changed  for  runs  4 and  5 where  waves  were 

introduced  simultaneous ly  with  the  start  of  the  tide,  i.e.,  part  A was 
eliminated. 
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t:  6 minutes 


t:  3 minutes 


high  water  slack 


t:  9 minutes  t:  12  minutes  t:  15  minutes 

Figure  16.  Time-lapse  photos  showing  evolution  of  a long  channel  bed 
(run  9A,  cycle  9)  under  tidal  action  (ebbtide).  Tidal 
period  is  30  minutes.  Numbers  are:  (1)  2- foot  barrier 
crest;  (2)  ocean:  (3)  waterline. 
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Period  and  height  (flap  stroke)  were  constantly  checked  at  regular 
intervals  throughout  the  run  to  ensure  uniform  wave  conditions.  Wave 
heights  varied  slightly  because  of  the  varying  water  depth  (submergence 
of  flap)  and  the  presence  of  currents  from  or  into  the  inlet.  Mild  waves 
(deepwater  steepnesses  smaller  than  0.0215)  were  selected  for  most  runs; 
only  a few  steep  waves  were  studied  because  of  time  limitations. 

The  evolution  of  the  beach  and  inlet  channel  was  different  from  that 
observed  under  tidal  action.  The  beach  geometry  adjacent  to  the  inlet 
was  shaped  in  response  to  the  waves  and  currents,  the  sand  transport 
towards  the  bay  increased,  the  channel  banks  became  flatter,  and  ebbtide 
channels  became  more  defined.  The  inlet  region  (short  channels)  pre- 
sented the  funnellike  shape  commonly  observed  in  natural  inlets,  with  the 
narrowest  section  (throat)  on  the  bay  side.  This  region  initially  ad- 
vanced towards  the  bay  at  a fast  rate,  until  it  reached  a position  when 
the  bay  tide  became  periodic;  some  long  channels  tended  to  meander. 

A series  of  photos  in  Figure  17  show  the  inlet  configuration  and  the 
current  and  wave  patterns  for  a run  with  a short  channel;  note  the  sand- 
spits,  the  geometry  of  the  inlet  region, and  the  wave  action  in  the  bay 
during  floodtide. 

Water  levels  were  usually  recorded  continuously  from  the  start  to  end 
of  a run.  Current  velocities  at  the  inlet  and  extreme  water  levels  in 
the  bay  were  recorded  during  the  last  few  cycles.  The  carriage  support- 
ing the  instruments  was  positioned  so  that  the  gages  and  the  current 
meter  were  at  the  throat  region  at  all  times.  The  run  was  usually  com- 
pleted on  the  same  day;  however,  in  a few  instances  the  run  was  halted 
and  resumed  the  following  day.  Halting  the  run  was  necessary  for  long 
inlet  channels  (wide  barrier)  that  required  long  run  times.  The  dura- 
tion of  runs  nB  ranged  from  4 to  25  hours. 

As  in  part  A,  most  runs  ended  at  slack  water  following  the  ebbtide 
flow,  and  the  final  procedures  discussed  previously  were  repeated.  Inlet 
cross  sections  were  spaced  at  0.25-  to  0.20-foot  intervals. 

After  the  final  measurements  were  taken,  some  of  runs  B were  extended 
from  2 to  18  hours  and  subjected  to  the  same  environment  (run  C)  to  detect 
further  evolution  and  movement  of  the  throat  region  (run  8C) , hydraulic 
changes  due  to  steeper  waves  (run  11C),  or  both  the  effect  of  time  and  of 
ending  the  run  after  the  floodtide  flow  on  the  cross-sectional  area  of 
the  throat  (run  12C).  Two  exceptions  to  the  standard  procedure  were: 

(a)  Run  18B,  for  which  the  channel  from  run  17B  was  used  as  a pilot 
channel,  and  (b)  tidal  period  increased  from  30  to  45  minutes  and  waves 
were  introduced  immediately  after  the  start  of  the  run.  Run  9B  followed 
run  9A  without  interruption. 

f . Inlet  Controlled  by  Jetties  (Runs  7C  and  18C) . Tests  with 
jetties  followed  runs  7B  and  18B  with  the  wave  conditions  and  ocean  tide 
unchanged.  Idealized  jetties  were  formed  using  two  long  thin  aluminum 
sheets  riveted  to  a reinforced  frame  (Fig.  18)  which  was  installed  in 
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the  channel  from  the  previous  run.  The  vertical  walls  were  parallel 
in  the  straight  section  and  alined  to  the  initial  pilot  channel.  The 
jetties  extended  into  the  ocean  offshore  of  the  low  water  level  mark  on 
the  beach.  The  bed  between  the  walls  was  undisturbed.  Lengths  of  the 
jetty  straight  sections  were  12  and  6 feet,  and  width  of  the  channels 
2 and  1.6  feet,  respectively. 

Before  the  start  of  the  run  the  superelevation  measured  in  part  B was 
introduced  in  the  bay,  separated  from  the  ocean  by  a semicircular  barrier 
placed  on  the  bay  beach.  Waves  were  introduced  immediately  after  the 
floodtide  flow  started. 

The  first  two  photos  in  Figure  18  show  the  "uncontrolled"  inlet  bed 
(run  18B)  before  placement  of  the  jetties  and  the  channel  bed  at  the  end 
of  the  "controlled"  inlet  run.  The  photos  in  this  figure  illustrate  the 
flow  and  wave  patterns  at  the  inlet  throughout  a cycle,  e.g.,  the  large 
amount  of  energy  transmitted  into  the  bay  by  the  waves  traveling  along 
the  channel  (smooth  walls)  during  floodtide. 

Durations  of  these  two  runs  were  24  hours  for  run  7C,  and  7 hours 
for  run  18C. 

g.  Freshwater  Discharge  into  Bay  (Run  16C) . After  cross-sectioning 
and  profiling  of  the  channel  for  run  16B,  the  test  was  resumed  under  the 
same  wave  and  tide  conditions.  One  and  one-half  cycles  later  a constant 
"freshwater"  discharge,  from  the  city  water  supply,  was  introduced  into 
the  bay  (at  the  comer  farthest  from  the  inlet)  until  the  end  of  the  run. 
This  discharge  was  measured  by  a rotameter.  The  constant  rate  of  flow 
was  10.2  gallons  per  minute  which  corresponded  to  about  9 percent  of  the 
maximum  discharge  through  the  inlet  at  strength  of  ebbtide  for  run  16B. 

A comparison  of  the  configuration  of  the  channel  bed  before  the 
start  and  after  the  end  of  this  run  which  lasted  4.5  hours,  is  shown  in 
Figure  19.  The  main  difference  is  the  shape  and  size  of  the  ripples 
close  to  the  throat,  probably  caused  by  the  increase  in  duration  of  the 
ebbtide  flow  and  the  superelevation  and  the  peak  velocity  at  the  throat. 
The  minimum  cross-sectional  area  below  MSL  was  equal  to  that  measured 
before  the  start  of  the  run. 

h.  Fixed  Ripple  Bed  under  Tidal  Action  (Runs  12F.  and  12F).  Two 
similar  runs  (12A  and  -12A)  were  conducted  as  described  for  runs  nA  to 
determine  the  variation  of  the  friction  characteristics  of  a short  inlet 
resulting  from  the  orientation  of  the  bed  ripples.  The  resulting  channel 
beds  were  fixed  after  the  end  of  ebbtide  and  floodtide  flows.  Tidal 
action  was  resumed  afterwards.  Fixing  of  the  channel  bed  was  done  by 
sprinkling  fast-setting  plaster  on  the  wet  ripples  and  beach;  thus,  a 
0.004-foot-thick  fixed  plaster  layer  covered  the  entire  inlet  region. 
Figure  20  shows  the  two  fixed  channels;  run  12E  is  on  the  ebbtide  ripple 
bed  and  run  12F  is  on  the  floodtide  ripple  bed.  The  pilot  channel  trace 
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Bed  after  ebbtide  (Ran  12K) 


Bed  after  floodtide  (Run  12F) 


figure  20.  I’hotos  showing  short  and  deep  inlet  beds  "fixed"  with 
plaster  after  last  ebbtide  and  after  last  floodtide, 
under  tidal  action  only.  Numbers  are:  (1)  sand  barrier 

crest;  (2)  distance  equals  1 foot  for  photo  scale; 

(3)  ocean  section  of  pilot  channel  unaffected  by  currents 

(4)  edge  of  plaster  layer. 


on  the  oceanside  is  still  visible  in  the  figure  because  the  cross- 
sectional  area  of  the  channel  is  large  (channel  depth  large  compared  to 
tidal  range),  and  the  currents  are  not  strong  enough  to  move  sediment  in 
this  region.  The  plaster  layer  thickness  was  considered  in  the  compu- 
tation of  areas,  widths,  and  perimeters. 

Duration  of  each  run  was  4 hours.  The  plaster  layer  in  run  12F  was 
carefully  removed,  and  the  initial  movable  bed  preserved  as  the  pilot 
channel  for  run  12B  with  wave  action. 

i.  Fixed-Bed  Prismatic  Channel  (Run  11).  A special  test  was  con- 
ducted on  a trapezoidal  inlet  (1.70- foot  base;  30°  banks)  cut  in  a nar- 
row sand  barrier  (bottom  length  7 feet)  and  fixed  with  plaster  to  con- 
trast the  results  from  the  movable-bed  studies  with  those  of  idealized 
fixed  prismatic  channels. 

The  test  involved  only  tidal  action  (run  11A)  during  a 3-hour  period, 
followed  by  4 hours  of  simultaneous  tidal  and  low  wave  action  (run  1 1 B) . 
Run  11C  used  tidal  and  steep  wave  actions  for  2 hours.  In  each  run,  the 
time  required  to  establish  a periodic  bay  tide  was  considerably  shorter 
with  the  fixed-bed  model  than  the  previous  movable-bed  tests. 


EXPERIMENTAL  DATA 


1.  Measurements. 


a.  Ocean  and  Bay  Tides  and  Current  Velocities.  Samples  of  the 
charts  used  to  record  tides  in  the  ocean,  inlet  channel  and  bay  (resist- 
ance gages) , and  flow  velocities  (current  meter  at  about  60  percent  of 
the  depth)  at  the  throat  or  at  midchannel  for  runs  A,  are  shown  in  Figure 
21.  Tide  scales  were  not  necessarily  the  same  for  the  three  tide  gages-- 
the  signal  from  the  bay  was  amplified  two  times  for  better  interpretation, 
especially  in  the  runs  with  small  bay  ranges.  The  figure  shows  slack 
watermarks  at  the  bottom  of  the  chart.  Chart  speed  was  5 millimeters  per 
minute  and  the  galvanometer  error  +0.25  millimeters  over  the  central  4 
centimeters . 


In  runs  B and  C,  the  wave  action  at  the  throat  was  sometimes  so  high 
during  floodtide  that  the  water  level  gage  recorded  the  wave  spectrum. 
These  high-frequency  signals  were  eliminated  for  the  remaining  tests  by 
installing  filters  in  the  electronic  circuits  (300-microfarad  capacitors) 
that  allowed  recording  of  the  mean  water  elevations  during  the  run. 

Ocean  and  bay  tides  were  also  recorded  by  separate  float  recorders. 

A pen  was  mounted  on  a rod  attached  to  a float  to  give  a scale  of  1:1 
and  a chart  speed  of  15  millimeters  per  minute.  A comparison  of  readings 
made  with  the  micro manometers  and  the  float  recorders  showed  a difference 
of  less  than  4 percent  of  the  average  reading. 

b.  Water  Surface  Slope  at  the  Inlet.  The  remaining  tests  included 
recording  of  the  water  elevations  and  determination  of  the  water  slope 
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Figure  21.  Samples  of  current  velocity  and  ocean,  inlet,  and  bay 
tide  recordings  by  resistance  gages.  Last  cycles  of 
tidal  action  (run  15A)  and  tidal  and  wave  actions 
(run  15B). 
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in  the  inlet  throat  region.  Typical  records  are  shown  in  Figure  22; 
traces  in  the  lower  part  of  the  figure  (run  15B)  are  the  water  surface 
difference  in  elevation  between  the  locations  of  the  two  resistance  gages 
(separated  3 to  4 feet  along  the  channel  centerline) , and  the  piezometric 
difference  of  the  same  points  obtained  by  placing  normal  resistance  gages 
each  inside  two  1-inch  diameter  stilling  wells  connected  by  thin  plastic 
tubing  to  intakes  placed  in  the  sand  about  3 inches  above  the  basin 
bottom  and  directly  below  the  resistance  gages.  The  stilling  wells  were 
4 inches  apart  and  always  surrounded  by  bay  water  so  that  any  change  in 
water  temperature  was  the  same  for  both,  and  the  signal  was  not  affected 
by  temperature  change.  Water  stage  at  the  throat  (located  midway  be- 
tween the  slope-measuring  gages)  was  also  recorded  on  the  same  chart. 

The  upper  part  of  Figure  22  shows  the  traces  corresponding  to  the  tide- 
only  condition. 

c.  Surface  Waves.  The  two  wave  gages  were  installed  along  the 

beach  centerline,  one  at  the  beach  toe  and  the  other  3 to  4 feet  ocean- 

ward.  Waves  were  modified  by  the  currents  and  water  depths  (tide) , and 

by  reflections  and  uncontrolled  vibrations  from  the  wave  machine.  Chart 
speeds  were  5 and  0.5  millimeters  per  second. 

d.  Flow  Cross  Sections.  Typical  sections  for  a short  and  deep 
channel,  and  location  of  the  current  meter  propeller,  are  shown  in 
Figure  23.  The  trace  of  the  channel  bed  surface  at  different  stations 
was  obtained  from  the  planigraph  used  to  reproduce  the  cross  sections. 

e.  Overhead  Photography.  Vertical  photography  was  taken  to  docu- 

ment all  inlet  runs.  Photos  in  Figure  24  show  the  channel  bed  (long 
channel)  at  the  end  of  a run  under  tidal  action  only,  at  the  end  of  the 
run  that  followed  (tidal  and  wave  action)  and  at  the  end  of  a run  exten- 
sion. Runs  shown  in  the  figure  ended  at  slack  water  after  ebbtide  flow, 

as  indicated  by  orientation  of  the  sand  ripples.  In  a few  cases,  the 
channel  was  contoured. 

f.  Inlet  Channel  Profile.  The  profile  along  the  centerline  of  a 
pilot  channel  (Fig.  25)  was  obtained  by  computing  depths  below  MSL  from 
readings  (to  the  thousandth  of  a foot)  of  the  bed  taken  with  a point  gage 
mounted  on  the  carriage.  It  covered  the  distance  from  the  bay  delta  to 
the  small  ripples  found  on  the  original  beach.  Care  was  taken  to  record 
the  lowest  and  highest  points 'of  the  channel  ripples  as  well  as  their 
shapes . 

g.  Logbook.  High  and  low  water  and  SWL  readings  from  the  micro- 
manometers, water  temperatures,  hook  gage  readings,  inflow  manometer 
readings,  current  meter  locations,  and  all  other  pertinent  information 
(including  visual  observations)  were  registered  in  a logbook  during  the 
experiments . 

2.  Measurement  Problems. 

Temperature  variations  between  the  bay  and  ocean  occasionally  caused 
significant  jumps  in  electronically  measured  water  levels;  these  records 
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Sample  recordings  of  water  surface  slopes  at  midchannel 
(run  15A)  and  at  the  throat  (run  15B) . Curve  at  center 
is  proportional  to  the  slope  given  by  piezometers. 


section  for  a short  and  deep  channel. 


channel.  Numbers  are:  (1)  low  waterline; 

(4)  0.04  foot  below  SWL. 


were  corrected  by  manually  measuring  the  highs  and  lows.  Therefore,  some 
slope  measurements  and  estimates  of  timelags  between  high  and  low  waters 
are  only  approximate  (run  15B) . 

Current  velocity  measurements  may  also  contain  slight  errors  due  to 
complex  flow  induced  by  bed  forms  and  the  variable  depth  of  the  current 
meter  throughout  the  tidal  cycle. 

3.  Reduction  and  Presentation  of  Data. 

Since  drift  of  the  resistance  gage  signals  was  uncontrollable,  the 
resulting  records  were  used  only  for  the  measurement  of  time  durations 
and  lags.  Float  gages  were  used  to  determine  MSL  and  MBL  (mean  bay 
level),  as  well  as  other  parameters  involving  water  elevations. 

Figures  26  and  27,  and  the  list  of  Symbols  and  Definitions  assist 
in  understanding  the  meaning  of  the  various  parameters  given  in  the 
tabulation  of  data. 

Elevations  were  read  from  the  float  records  by  using  a special 
eyepiece  that  permitted  readings  to  0.0005  foot.  The  elevations  from 
the  resistance  gages  (at  the  inlet  channel)  were  made  by  using  the 
appropriate  calibrations.  Tidal  ranges  in  the  bay  and  ocean  were  com- 
pared to  those  given  by  the  micromanometers;  the  micromanometers  were 
taken  as  the  most  accurate,  only  for  very  small  bay  ranges.  All  of 
these  readings  corresponded  to  the  last  cycle  of  the  run. 

Lags  , duration,  and  the  measurement  of  slope  of  the  bay  tide  curve 
at  maximum  tlow  (d’n^/dt)  were  taken  from  the  resistance  gage  charts. 

The  slope  was  often  compared  to  measurements  from  the  float  charts  and, 
at  times,  the  charts  were  chosen  as  definitive  for  computations  because 
of  irregularities  or  difficulty  in  reading  the  bay  tide  curve.  Justi- 
fication for  compilation  of  data  at  maximum  flow  is  given  in  Section  V. 

The  flow  cross-sectional  areas  were  obtained  by  planimetering.  The 
Ac  > WMSX,  and  p values  were  selected  as  follows.  Since  the  shapes 
and  areas  of  the  cross  sections  seldom  change  from  station  to  station, 
average  values  were  computed  for  tidal  action  only  (runs  A).  Whereas, 
for  runs  with  waves  (a  well  defined  cross  section  with  a minimum  area) 
the  throat  characteristics  were  used  as  data  for  the  computations  in 
Section  VI. 

Data  on  the  geometric  characteristics  of  the  cross  sections  are 
given  in  normalized  form  in  Section  V and  in  the  Appendix. 

Wave  heights  measured  at  the  toe  of  the  beach  were  averaged  over 
several  tidal  cycles  close  to  the  end  of  the  run  and  then  transferred  to 
deep  water  (linear  theory)  to  obtain  wave  height,  length,  steepness,  and 
wave  power. 

Hydraulic  and  other  pertinent  data  are  presented  in  Tables  1 to  4, 
and  includes  the  main  variables  in  the  ocean,  inlet  channel,  and  bay. 


62 


Channel < , Ocean > . Throat 


Elevation  Negative 
Below  MSL 


Figure  26.  Schematic  representation  of  simultaneous  water  elevations  in 
the  ocean,  inlet,  and  bay;  lags  and  durations;  and  the 
velocity  at  the  inlet  throat  over  a tidal  cycle. 
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b)  A-A  Section:  Typicol  Flow  Cross  Section  Perpendiculor  to 
to  Channel  Alinement  ( overage  for  runs  A and  at  throat 
for  runs  8 and  C ) 


Figure  27.  Diagrams  identifying  the  main  parameters  of  the  sand 
barrier  profiles  and  the  flow  cross  sections. 
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V.  CORRELATION  OF  HYDRAULIC  VARIABLES 


Using  the  data  in  Tables  1,  2,  and  3,  each  experimental  run  was 
identified  by  its  characteristic  coefficients  of  repletion  (different 
for  ebbtide  and  floodtide)  as  defined  by  Keulegan  (1967,  equation  4)  to 
obtain  a direct  comparison  with  his  theory.  Corresponding  values  of  the 
energy  losses  coefficient,  F,  were  also  computed. 

Subsequently,  and  based  on  measurements,  computations  were  made  of 
the  tidal-damping  ratio  (a^/ap);  superelevation  (A/Rc)  where  A is 
mean  bay  level  above  MSL;  timelag  between  high  waters  (e^) ; dimen- 
sionless maximum  mean  flow  velocities  (V^,) ; absolute  values  of  tidal 
prisms,  ft;  inlet  channel  geometric  and  hydraulic  characteristics 
(Ac,  R-MSL>  D#) ; maximum  mean  flow  velocities  at  ebbtide  and  floodtide 
(V maxe>  vmaxf)>  and  maximum  discharge  (Q^^) . 

The  results  are  given  in  Tables  4,  5,  and  6. 


Computations  of  Parameters  from  Experimental  Data. 


a.  Computations  of  K.  The  Keulegan  (1967)  assumptions  of  a vertical 
bay  wall  and  uniform  bay  level  motion  were  observed  to  hold  for  the  inlet 
conditions  tested;  therefore,  the  Keulegan  repletion  coefficient  is  com- 
puted using  ebbtide  and  floodtide  conditions  separately  to  determine 
Kg  and  Ky\ 


Since  the  bay  basin  walls  are  vertical  (6=0),  the  instantaneous 
rate  of  filling  (floodtide  flow)  or  emptying  (ebbtide  flow)  of  the  bay 
is  given  by: 


Q = A b 


d hb 

dt  ’ 


(12) 


dh^/dt  being  the  slope  of  the  bay  tide  curve  (values  in  Table  3 are  the 
maximum  values  of  these  slopes).  This  rate  is  equal  to  the  flow  dis- 
charge through  the  inlet  channel  (conservation  of  mass). 

At  maximum  discharge  the  acceleration  term  8V/8t  in  equation  (7) 
is  assumed  to  be  negligible.  Therefore,  equation  (3)  can  be  used  to 
compute  the  corresponding  values  of  K,  i.e., 
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for  floodtide,  and 


Kg  = 


/,2Rr  /(h&  - h0)( 


A^ 


dh& 


jnax 


dt 


(14) 


for  ebbtide,  where  subscript  f is  for  floodtide  and  e is  for  ebb- 
tide. ka  is  the  dimensionless  cross-sectional  area  (Ae/Ae  at  MSL)  used 
in  these  two  equations  to  account  for  the  varying  cross-sectional  area 
with  tide  stage. 


Values  of  Kf  and  Kg  for  each  run  were  computed  with  equations  (13) 
and  (14),  by  using  all  the  respective  measured  parameters  T,  Rg,  h£>, 
h o,  [dhfc/dt]  max  (from  bay  tide  curves)  and  where  the  ratio  of  the 
actual  channel  cross-sectional  area  at  maximum  discharge  (e.g.,  at  max- 
imum floodtide  discharge)  to  the  area  at  MSL  is: 


. i * (WA 1SL  + Z hif)  hif 

A of  = — , (15) 

Ao 

with  A a the  average  area  at  MSL,  ViMSL  the  width  at  MSL,  £ the  bank  slope 
(horizontal  distance  per  unit  of  vertical  distance)  and  h if  the  inlet 
water  level  above  MSL  at  maximum  floodtide  discharge. 


b.  Computation  of  F.  The  coefficient,  F,  accounts  for  the  unit 
energy  losses  at  maximum  inlet  discharge  due  to  friction  (fL/4R) , flow 
contraction  into  the  inlet  channel  (entrance  losses,  k en) , flow  expan- 
sion into  the  bay  or  ocean  (exit  losses,  kex) , and  other  minor  losses 
(included  in  ken  + kex) . Values  of  F may  be  computed  both  at  max- 
imum ebbtide  and  floodtide  discharges  using  Keulegan's  definition  of  K: 
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c.  Computation  of  k/f.  The  product  K/F  is  a dimensionless 
number,  constant  throughout  a tidal  cycle  for  a given  inlet,  where: 


K/F  = Kf/rf  = Ke/F7  = JEL  . 

11  e A, 


(18) 
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K/f  is  a useful  variable  since  its  computation  is  based  only  on  non- 
friction terms,  measurable  both  in  the  laboratory  and  in  the  field.  For 
this  parameter  it  is  unnecessary  to  estimate  the  friction  coefficient,  f, 
the  exit  and  entrance  energy  loss  coefficients  and  the  inlet  length,  and 
the  effect  of  the  ocean  tide  on  the  channel  is  reflected  in  the  value  of 
the  inlet  cross-sectional  area  for  an  inlet  that  has  reached  an  "equilib- 
rium" (for  these  tests,  values  of  K^F  were  computed  when  the  bay  tide 
became  periodic). 

d.  Computation  of  Dimensionless  Current  Velocity  and  Discharge. 

The  maximum  discharge  through  the  inlet  channel  is  given  by: 


The  maximum  slope  of  the  bay  tide  lasts  long  enough  to  validate  the 
assumption  that,  for  all  practical  purposes,  the  maximum  mean  velocity 
across  a section  is  given  by: 


v = ^rncx 
max  A a 


From  O'Brien  and  Dean  (1972),  dimensionless  values  of  maximum  mean 
current  velocities  at  the  selected  cross  section  were  determined  from 
experimental  results  for  both  ebbtide  and  floodtide  using  the  defini- 
tion: 


Qrnax^ 
ttA b Ro  ’ 


or  using  equation  (18): 


v'  = V a* 

mn.T.  v ivn'r  r ' ' 


Dimensionless  values  of  maximum  discharge  were  also  computed  from 
the  definition: 


^max  = CD  • Ae  • /gR0 


which  yields,  combined  with  equation  (20),  the  "discharge  coefficient,' 
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and,  using  equation  (22), 


1 " 


CD 


\max 

K/p  a;  ‘ 


(25) 


The  absolute  theoretical  velocity  is,  from  equation  (24), 


max 


(26) 

i 


which  may  differ  from  the  experimental  values  computed  in  equation  (20). 


e . Computation  of  Friction  Factor  (f)  and  Roughness  Coefficient  (n) . 
Estimations  of  friction  factors  at  maximum  discharge  were  made  using  the 
channel  length,  L,  and  the  cross  section  as  shown  by  the  measurements  in 
Table  2,  with  the  friction  factors  for  ebbtide  and  floodtide  flow  defined 
by: 


f L (F  " {^en  * ^ex))  . (27) 

For  these  calculations,  it  is  assumed  that  + k ex  = 1-0  (Huval 
and  Wintergerst,  in  preparation,  1977). 


Since  the  relationship  between  f and  the  Manning  roughness 
coefficient,  n,  is  given  by  equation  (11),  estimates  of  n may  be 
obtained  by: 

n = 0. 18S  • R2^  • /|p-l) / L . (28) 

Experimental  values  of  f range  from  about  0.01  to  0.65,  but  most 
range  from  0.05  to  0.20  (Tables  4,  5,  and  6). 

Computed  values  of  n range  from  about  0.006  to  0.049,  but  most 
range  from  0.013  to  0.030. 

3.  Computation  of  Theoretical  Parameters  by  the  Lumped  Parameter 
Approach. 

The  computer  program  in  Huval  and  Wintergerst  (in  preparation, 

1977)  was  modified  to  use  the  experimental  data  collected  in  this 
investigation.  The  input  to  the  new  program  (DATA THE)  included: 

(a)  All  measured  data  necessary  to  compute  K and  F at  the 
maximum  ebbtide  and  floodtide  discharge; 


(b)  a sinusoidal  ocean  tide  with  range,  R0,  and  period,  T, 

(c)  initial  values  of  man  flow  velocity  and  bay  water  ele- 
vation. 
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(d)  coefficients  indicating  the  time  interval  selected  for 
integration  of  the  differential  equations  and  the  number  of  tidal 
cycles  necessary  for  the  solution  to  reach  a periodic  bay  tide, 
and 

(e)  3=0  (vertical  bay  walls). 

The  program  determines  n for  ebbtide  and  floodtide  flows;  then, 
two  different  integrations  are  made,  one  for  each  of  the  n values.  It 
is  assumed  that  the  roughness  coefficient  at  maximum  ebbtide  or  floodtide 
discharge  is  unchanged  throughout  the  process  of  integration.  Therefore, 
for  each  run  two  similar  prismatic  channels  of  length  and  cross-section 
characteristics  are  used  (Table  2),  but  with  different  roughnesses. 

The  integration  procedure  is  explained  by  the  description  of  program 
INLET  (Huval  and  Wintergerst,  in  preparation,  1977).  Program  output  in- 
cludes the  bay  tide,  the  mean  current  velocity,  the  discharge  at  the  in- 
let as  a function  of  time,  the  net  discharge  for  a tidal  cycle,  maximum 
and  minimum  water  levels,  timelags  and  durations,  elevations,  etc. 

(Table  7).  These  results  constitute  the  theoretical  values  (lumped 
parameter  approach)  that  will  be  compared  to  the  experimental  values  in 
Tables  4,  5,  and  6.  Figures  at  the  bottom  of  Table  7 correspond  to  theo- 
retical values  of  the  parameters  used  for  comparison  to  experimental 
results.  Therefore,  two  sets  of  solutions,  one  for  each  of  the  two 
roughness  coefficients  computed,  are  obtained. 

After  a thorough  study  of  the  results,  it  was  decided  that  a repre- 
sentative lumped  parameter  approach  would  be  the  solution  for  the  ebb- 
tide flow  roughness  coefficient,  which  implies  that  theoretical  results 
correspond  to  values  of  F and  K at  maximum  ebbtide  flow  only.  The 
K values  covered  a range  wider  than  that  corresponding  to  the  floodtide 
flow  and  yielded  more  consistent  results. 

Computer  programs  used  in  the  calculations  are  listed  in  Mayor-Mora 
(1973). 

4 . Comparison  of  Measured  and  Theoretical  Tides  and  Currents. 

The  theoretical  values  obtained  by  the  lumped  parameter  approach  of 
ocean  and  bay  elevations,  the  mean  current  velocity  at  the  inlet,  and 
the  rate  of  discharge  using  the  coefficient  of  repletion  defined  by 
ebbtide  flow  (run  8A)  are  shown  in  Figure  28.  The  difference  in  bay 
level  when  K is  defined  by  floodtide  is  shown  in  Figure  29  (run  8A) ; 
the  figure  includes  theoretical  bay  levels  defined  by  ebbtide  flow  when 
waves  and  tides  are  present  (run  8B) . 

A comparison  of  measured  data  and  lumped  parameter  theoretical  results 
gave  the  best  agreement  for  a coefficient  of  repletion  defined  for 
ebotide  flow.  Kg;  therefore,  this  parameter  was  used  in  analysis. 
Figures  30  and  31  show  the  agreement  between  measured  and  predicted 
hydra-.il ic  parameters  for  ebbtide  flow  prediction  with  tides  only  and  for 
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Run  Run 

8B  0.397  0.088  0.223  317.8  329  4 0 276  0 1255  0 1333  0.567  0.176  0 243  8B  0 202  0122  0 605  161  4 329  4 0 276  0 1255  0.1333  0.581  0.192  0 273 

8B  0.331  1.55  0.412  1 63  0 0078  0.2877  0.402  0.014  0 0061  0 699  0.302  0 98  88  0.158  6 83  0 412  0.77  0 0250  0 1353  0.189  0 145  0 0061  0 329  0 312  0 95 

8B  0.433  0.331  1.55  0.412  1.55  0 0078  0.4109  0.574  0.014  0 998  0.98  8 B 0.419  0.158  6.83  0.412  0.80  0 0250  0.2152  0 301  0.145  0.523  0.95 
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action  only  and  ebbtide  flow  coefficient  of  repletion,  K=0.179). 


OZOO 
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as  functions  of  time.  Comparison  of  measured  values  and  theoretical 
solution  by  the  lumped  parameter  approach;  low  ebbtide  repletion 
coefficient,  K=0.158  (with  waves). 


tide  and  wave  conditions  for  small  values  of  K„  (runs  8A  and  8B,  respec- 
tively) . Experimental  and  theoretical  results  for  relatively  large  values 
of  Kg  are  compared  in  Figures  32  and  33  (runs  12A  and  12B,  respectively) 

A possible  reason  for  the  better  agreement  between  observed  and  pre- 
dicted parameters  for  K defined  at  ebbtide  is  that  during  ebbtide  the 

Keulegan  assumption  of  open  channel  flow  is  better  approximated  than 
during  floodtide  when  the  inlet  acts  more  as  an  orifice. 

It  was  also  found  that  wave  action  in  addition  to  tides  reduces 

agreement  between  predicted  and  measured  water  levels  and  velocities 
(Figs.  30  to  33).  This  may  be  due  to  wave-induced  velocities  or  a 
change  in  the  effective  geometric  characteristics  of  the  inlet. 

5.  Dimensionless  Parameters  as  Functions  of  Kg. 

Five  dimensionless  inlet  hydraulic  parameters  are  presented  as  func- 
tions of  inlet  repletion  coefficients  (Figs.  34  to  39).  The  parameters 
are  plotted  using  open  symbols  for  tide  only  experimental  tests,  and 
closed  symbols  for  experimental  results  from  tests  with  tides  and  waves. 
Keulegan's  (1967)  results  are  shown  by  a solid  line,  and  results  from 
the  lumped  parameter  techniques  (Huval  and  Wintergerst,  in  preparation, 
1977)  by  a dashline. 

The  tidal  damping  coefficient,  ab/ao,  is  plotted  as  a function  of 
Keulegan's  repletion  coefficient  for  ebbtide,  Kg,  in  both  experimental 
results  and  from  theory  (Fig.  34).  All  results  are  similar  and  lead  to 
the  following  conclusions  for  the  parameter  a£,/a 0'  (a)  Inertia  is  not 
important  in  bay  tide  levels  for  these  conditions  because  Keulegan  (no 
inertia)  and  lumped  parameters  (inertia)  techniques  give  similar  results, 
(b)  ab/ao  is  the  same  for  tests  with  and  without  waves,  and  (c)  there  is 
excellent  agreement  between  theoretical  and  experimental  results. 

An  alternate  method  of  presenting  the  tidal  damping  coefficient  is 
to  plot  the  variable  against  the  inlet  parameter,  Kgv'F  (Fig.  35).  The 
same  trends  and  good  agreement  are  found  between  ab/ao  and  Kg^F  as  with 
the  repletion  coefficient;  therefore,  either  inlet  parameter  may  be  used. 
Since  this  was  found  for  other  dimensionless  hydraulic  variables  as  well, 
only  plots  of  Ke  are  presented. 

Dimensionless  lags  between  high  and  low  waters,  equal  to  the  timelag 
divided  by  the  tidal  period,  show  deviations  from  the  Keulegan  predic- 
tions and  best  agreement  with  the  lumped  parameter  solution  (Figs.  36 
and  37) . Experiments  show  that  lags  between  high  water  are  shorter  and 
lags  between  low  water  are  longer  than  Keulegan's  predictions.  This  may 
be  due  to  the  change  in  the  cross-sectional  area  and  the  size  of  the 
friction  and  inertia  terms  over  a tidal  cycle. 

The  dimensionless  maximum  inlet  velocity,  Vrnax>  from  experiments 
shows  good  agreement  with  Keulegan's  theory  (Fig.  38). 
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Keulegan's  repletion  coefficient  for  ebbtide. 


Dimensionless  timelag  between  high  waters  as  a function  of  Kg- 


SCO 


Figure  37.  Dimensionless  timelag  between  low  waters  as  a function  of  K( 


Dimensionless  maximum  inlet  velocity  as  a function  of  Kg. 


Figure  39.  Dimensionless  bay  superelevation  as  a function  of  K 


Dimensionless  bay  superelevation,  A/ ao,  (superelevation  divided  by 
ocean  tide  amplitude)  is  plotted  in  Figure  39.  Keulegan  predicts  no 
superelevation  which  the  tests  show  for  inlets  with  Ke  greater  than  1. 
However,  if  Ke  is  less  than  1,  bay  superelevations  may  be  significant 
with  superelevation  higher  in  tests  with  waves  and  tides  than  in  tests 
with  only  tides.  High  bay  superelevations  may  be  partly  due  to  scale 
effects  of  the  relatively  large  model  wind  waves. 

6.  Effect  of  Ripple  Orientation  on  Friction  Coefficients. 


Above  a certain  critical  velocity,  inlet  bed  material  will  move  con- 
tinuously and  form  asymmetrical  ripples  with  the  lee  side  of  the  ripples 
generally  steeper  than  the  side  facing  the  current  (Figs.  40  and  41). 

To  study  the  effect  of  ripple  orientation,  tests  were  run  (Section 
III)  to  qualitatively  examine  changes  in  the  friction  factor,  f,  for 
different  flow  conditions.  Where  f is  defined  as: 


f = 4R 
L 


(F-l) ‘ 


with  ken  + kex  =1.0  assumed.  A periodic  tide  was  run  with  a short,  deep 
inlet  until  equilibrium  was  reached;  then  the  bottom  was  fixed  on  an 
ebbtide  plan  by  applying  a thin  layer  of  plaster.  Tests  were  continued 
until  f could  be  estimated  for  ebbtide  and  floodtide.  The  test  was 
then  repeated  by  fixing  the  inlet  with  ripples  formed  during  floodtide 
flow. 

The  geometric  properties  of  the  inlet  with  ebbtide-oriented  ripples 
are  the  same  as  those  with  floodtide  ripples  with  K/F  having  only  a 
4 percent  difference.  Since  the  bed  material  was  plaster  and  not  sand, 
these  tests  can  only  be  compared  with  each  other. 

The  tests  indicate  that  for  floodtide  ripples,  both  ebbtide  and 
floodtide  flows  at  maximum  discharge  have  the  same  friction  factor,  f; 
however,  for  ebbtide  ripples,  ebbtide  flow,  f,  was  higher  than  for 
floodtide  flow  (Fig.  42). 

This  change  in  friction  suggests  that  during  a tidal  cycle,  friction 
will  be  high  for  ebbtide  flow  on  ebbtide  ripples.  When  the  tide  changes 
and  ripples  have  not  had  a change  to  reorient,  floodtide  flow  will 
experience  low  friction  on  these  ebbtide  ripples.  This  difference  is 
partly  due  to  the  definition  of  friction  and  experimental  error,  and 
further  work  is  suggested. 

7.  Tidal  Inlet  Prism  versus  Area  Relationships. 

Tidal  prisms,  Q,  and  cross-sectional  areas,  Ac,  measured  experi- 
mentally in  these  tests  are  plotted  in  Figure  43  on  log-log  paper.  A 
best-fit  regression  analysis  shows  that  the  logarithms  of  prism  and 
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Typical  rippled  beds  after  ebbtide  flow  for  high  (F)  inlets 
(8)  Oceanside;  (1)  channel  centerline;  (2)  Stillwater  level 


Numbers  are 


0 Brien  (1969)  Ac  = 2.0X  I0~5  & 
Extrapolated  Field  Results  for 
Inlets  without  Jetties 


/ Jarrett  (1976)  Extrapolated  Field  Results 
/ 95 Pet  Confidence  Limits 


/ . v° 

/ • /•<*> 

Tide  and  Waves  //  / 
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o Tide  Only 
• Tide  and  Waves 


All  Tests 

A r = 6.0  X IO-3n0  735 
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Figure  43.  Model  inlet  tidal  prism  versus  area. 


area  are  linearly  related 
0.90  for  tests  with  tides 

Table  8.  Experimental 

with  correlation  coefficients  greater  than 
and  tests  with  both  waves  and  tides  (Table  8). 

relations  between  inlet  prisms  and  areas. 

Type  of  tests 

Relation 

Tide  only  (15  runs) 

n = 1400  A 1-73  or  A = 1.49  x 10"2  fiO-578 
c o 

Tide  and  waves  (20  runs) 

fl  = 1050  A 1<30  or  A = 4.76  x 10‘3  fi0-769 
c o 

All  tests  (35  runs) 

Q,  = 1050  A 1,36  or  A = 6.0  x 10'3  fi0-735 
a a 

Model  tests  show  that  for  a given  tidal  prism,  a model  inlet  will 
have  a smaller  cross-sectional  area  for  a test  with  waves  and  tides, 
than  an  inlet  with  tides  only.  This  suggests  that  waves  carry  sand 
to  the  inlet  and  reduce  the  inlet  area. 

Comparison  of  these  model  tests  with  results  from  models  by  McNair 
(in  preparation,  1977)  shows  that  the  two  sets  of  tests  had  similar  ranges 
of  experimental  conditions  and  results  of  tidal  area  versus  prism. 

Comparison  of  model  results  with  field  data  from  the  east,  west,  and 
gulf  coasts  (Jarrett,  1976)  shows  that  for  a given  tidal  prism,  model 
tests  predict  an  inlet  area  an  order  of  magnitude  higher  than  the  95- 
percent  confidence  limit  of  extrapolated  field  data  (Fig.  43).  This 
suggests  that  either  the  inlet  area  versus  prism  relationship  may  not  be 
an  exponential  relation  for  the  entire  range  of  conditions  from  model  to 
prototype  or  that  scale  effects  are  important.  A possible  scale  effect 
in  the  tests  is  the  use  of  typical  size  quartz  sand  which  does  not  model 
prototype  barrier  sediment  conditions. 

VI.  SUMMARY  AND  CONCLUSIONS 

A series  of  36  tests  were  conducted  in  an  idealized  movable-bed 
tidal  inlet  model.  Sinusoidal  tides  and  model  waves  were  run  contin- 
uously at  the  model  inlet  until  a periodic  tide  was  established  in  the 
bay.  Parameters  measured  included  ocean  and  bay  hydraulic  conditions, 
inlet  velocities,  and  the  inlet  geometry. 

Results  of  these  tests  show  that  two  techniques  accurately  predicted 
idealized  model  inlet  hydraulics.  The  Keulegan  (1967)  inlet  repletion 
coefficient  technique  predicted  the  ratio  of  the  bay  to  ocean  tide  am- 
plitudes for  a wide  range  of  model  conditions  when  the  repletion  coef- 
ficient was  defined  for  ebbtide  flow.  The  lumped  parameter  method 
(Huval  and  Wintergerst,  in  preparation,  1977)  which  extends  the  Keulegan 
method  by  considering  inertia  and  variable  inlet  geometry,  gave  a good 
prediction  of  maximum  inlet  velocity,  bay  superelevation,  and  the  time- 
lag  for  high  and  low  waters. 
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Additional  tests  on  a fixed-bed  inlet  with  ripples  showed  that  both 
ebbtide  and  floodtide  flows  had  the  same  friction  factor  at  maximum 
discharge  for  f loodtide-oriented  ripples,  but  higher  friction  over  ebb- 
tide ripples  for  ebbtide  flow. 

Analysis  of  model  tidal  prisms  and  inlet  areas  for  movable-bed  models 
shows  that  as  the  tidal  prism  increases  the  inlet  area  also  increases  for 
prototype  inlets  (O'Brien,  1969;  Jarrett,  1976).  However,  the  constants 
in  this  equation  are  not  the  same  for  model  and  prototype  inlets.  The 
reason  for  this  difference  between  model  and  prototype  constants  is  that 
forces  were  scaled,  while  sediment  was  not,  so  that  scale  effects  become 
important. 

Waves  superimposed  on  tides  reduced  the  inlet  cross-sectional  area 
for  a given  model  tidal  prism  with  a 40  percent  reduction  for  the  smaller 
areas  tested.  This  suggests  that  waves  cause  a net  transport  of  sediment 
to  the  inlet  zone  which  reduces  the  inlet  area. 
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APPENDIX 


MODEL  GEOMETRIC  DATA 


Plots  showing  model  data  at  mean  sea  level 
as  a function  of  stations  along  the  inlet. 
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